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ABSTRACT 
The article presents an effective numerical-analytical method for solving problems of hydro-viscose-elasticity and gives examples of 
calculating the natural oscillations of thin-walled cylindrical shells in contact with a viscous fluid by this method. The linear torsional 
vibrations of a cylindrical viscoelastic shell in contact with a viscous fluid are considered. The equation of motion of a viscoelastic 
shell is obtained on the basis of the moment theory of shells and satisfies the Kirchhoff Love and S.P. Timoshenko hypotheses. Fluid 
motion is described by the Navier Stokes equations. The problem is solved numerically by the methods of orthogonal sweep and 
Muller. The spectrum of complex Eigen frequencies and modes of oscillations is obtained depending on the wave number. 
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1. INTRODUCTION 


Despite a wide-ranging study of the interaction of elastic systems, there are still issues to be solved. The question of the effect of 
wave formation and compressibility of a liquid has not been fully investigated. Calculation methods are often cumbersome and 
require simplification. The interaction of the viscoelastic systems themselves in fluid flows has been little studied. The phenomenon 
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of the propagation of a wavelike motion of a fluid in elastic cylindrical shells attracted the attention of researchers [1-6]. In these 


works devoted to wave processes, the system “elastic cylindrical shell - ideal fluid” uses classical and refined shell equations, 
considers the effect of radial and longitudinal inertial forces, takes into account the average density of a fluid or gas flow. In [7-9], an 
analysis of the laws of the wave process in an elastic shell with a viscous fluid is performed within the framework of the model of 
linearized hydrodynamic equations of a viscous compressible fluid. Unlike the others, the system “cylindrical shell (elastic or 
viscoelastic) and liquid” (more perfect or viscous) is considered as a dissipative non-uniform mechanical system [10-12]. 


2. PROBLEM STATEMENT 


The natural oscillations of a cylindrical shell of radius infinite along the length of a deformable (viscoelastic) cylindrical shell are 
considered R; with constant thickness ho, shell density 9 , Eo- instant modulus of viscoelastic shell, Poisson's ratio V o, filled into a 


viscous liquid with a density QO, in equilibrium. Shell vibrations when exposed to internal pressure P(-P, »—P>+P3) described by 


equations [8, 12, 19]: 


4 oe I-vy J. \l-v, Jo°a 
La fun ente ye tN 5g i 


a 
E, ot a) 


Here E, - instant modulus, U = ii(u, Ug,U ). the vector of displacements of the points of the middle surface of the shell, 


& 


and for Kirchgor — Love shells, it has dimension equal to three (u, =U, Ug=Vv, U,= w), and for Tymoshenko-type 
shells, the dimension of the vector is five; here, in addition to the axial, circumferential, and normal displacements, the angles of 
rotation of the normal to the median surface in the axial and circumferential directions are added [12]; R;, (t _ T)-core relaxation. 


The amplitudes of oscillations are considered small, which allows one to write down the basic relations in the framework of the linear 
theory. The system of linearized equations of motion of a viscous barotropic fluid can be represented as [12] 
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Here in the equations (2) o = AG, oh) - velocity vector of fluid particles, and P - density and pressure 


* * * 
perturbations in the fluid, P) and ao — density and speed of sound in a fluid at rest; V , - kinematic and dynamic viscosity 
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* * 
coefficients; for the second viscosity coefficient A accepted ratio A =- =u; Piz> Pi» Pg - components of the stress tensor 


in a liquid. Equations (1) are respectively kinematic and dynamic boundary conditions, which, by virtue of the thinness of the shell, 
will be satisfied on the middle surface (r = R). Relations (1) and (2) represent a closed system of hydro-viscoelastic relations for a 
cylindrical shell containing a viscous compressible fluid. So for shells that obey the Kirchhoff-Love hypothesis, joint shells and liquids 
that are harmonic in the axial coordinate z and exponentially damped in time, or harmonic in time and damped in z are subject to 
study. Expanding equations (2) and (3) in coordinate form, it is easy to see that relations (2) - (3) break up into independent 
boundary value problems: 


Torsional vibrations 
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Let the wave process be periodic in z and attenuate in time, then the real wave number k is given, and the complex frequency is 
the desired eigenvalue. The solutions of boundary value problems (2) - (6) for the main unknowns satisfying the constraints imposed 
above do not depend on time and the z coordinate should be sought in the form [14] 


ice) 
Py, = ye. cos(mp)e"*", 


m=1 
a . + 
Pre = as cos(mp)e"* , 
m=l 
= re + 
Pro a yt, sin(mg)e*"™, (7) 


m=1 
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i => ii,,(U,, cos(m@),V,, sin(mg), W,, cos(m@) Je ; 


G= > (Fm COS), F,, sin(m@), F 


cos(m@) jen 


zm 


where 0,,7,,T 9, Ce a FD, F, - amplitude complex vector - function; k is the wave number; C is the phase 


velocity; w is the complex frequency; m is the circumferential wave number (the number of circumferential waves), taking the values 
m = 1, 2, 3... In the case of m = 0, axisymmetric oscillations. This approach will allow us to search for a solution to the problem for 
each fixed value of the circumferential wave number m independently. After performing in (5) the change of variables (7), the 
resolving relations describing the stationary torsional vibrations of the shell-liquid system are formulated as a spectral boundary 
problem for a system of two ordinary differential equations 
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We first investigate the oscillations of the fluid in the walls. Equations (9) can be converted to a single equation for displacement v 
d°9, d8 Oo 1 * 
4+ —* + (-k* +i =-—)9, =0; v al (10) 
dr rdr vor Po 


The solution of equation (10) is 
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Where Jiand N; Bessel and Neumann functions, respectively, of the first kind and first order; A: and Bi- arbitrary constants. Given 


the immobility of the shell, to obtain the dispersion equation 


J(Ryf-k +42) =0 (12) 


@, =—i(v"k? + B?) (13) 


from where 


Here through £ , marked by the roots of the Bessel function, referred to R. As can be seen from formulas (12), (13), the proper 
motions are always a periodic in time, and the nodal points are fixed (phase velocity C.=0), at that time, the steady-state motion is 
oscillatory, and the nodal points move with the speed Cy, monotonically increasing from zero to infinity with increasing or viscosity 
Vv" . These characteristic features of viscous medium motion will manifest themselves in the following more complex examples. We 


now consider relations (9) for the case of the internal location of the fluid. This problem can be solved in the same way, using special 
functions and we have the dispersion equation 


2 * J 
a + — =t(Z 02) _») 9 (14) 
a a phR J,@) 


Which was first obtained in the work of A. Guzya [7]. New notation introduced here. 
p= Qo ‘h= we R,,|/-k on se = a 
O. 


Shell shear wave velocity: Jo-zero order Bessel function. 

The direct solution of equation (14) faces certain difficulties caused by the necessity of calculating the Bessel function of the 
complex argument. Therefore, we study (14) using asymptotic representations of these functions for small and large arguments z. 
The smallness of z occurs at low-frequency oscillations. According to known expansions Jo and J; power series 


2 2 
Zz Zz Zz 
J, =l1-— =...3J,(z) =—(-~—+...); (15) 
pI i @= GIS +) 


a 2, 0 
Holding in the expansions (15) only the first terms, we get —k~ + =a = 0 
a 
Dispersion equation of torsional vibrations of a dry shell or filled with an ideal fluid, keeping in (14) the first two terms in each, 
we have the equation 
2 


— er iO Le -iZ)=0 16 
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The root of which, for example, in the case of steady-state oscillations is determined by the expression 
1/2 
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A physical interpretation of equation (14) is provided below. We now consider the situation when z is large enough, which 


corresponds to high-frequency oscillations and low viscosity. Substituting (1), and additionally assuming smallness Vv" compared to 


we obtain an approximate dispersion equation, which is also given in [7] 


2 * ~ 7 
[+i 

ee 1+,/— =0 (18) 
- oO PRiae 


Hence, when aspiring viscosity v to zero (as well as @ to infinity), we have a trivial result — —> Q, which was obtained at 


ke 1 


small @ from equation (16). Equation (18) is unacceptable at high viscosities. In this case, the phase velocity C increases unlimitedly 
with increasing w. The considered example testifies to the inconsistency of various asymptotic estimates in the region of average 
oscillation frequencies. Thus, when analyzing wave processes by asymptotic methods in the first approximation, it is not possible to 
establish the limits of applicability of the formulas obtained, as well as to estimate the error of calculations. In this paper, for solving 
spectral problems, we use direct numerical integration of resolving relations of the type (9) using the orthogonal sweep method in 
complex arithmetic. 
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Figure 1 Dependency Re @ from wave number k 


Such an approach makes it possible to avoid the above difficulties associated with the calculation of the Bessel functions of a 
complex argument. Another advantage is due to the specifics of the orthogonal sweep method, which, thanks to the ortho- 
normalization procedure, allows one to solve strongly rigid systems with a boundary layer. As a result of the conducted numerical 
study, it was found that the problem of natural oscillations (9) admits no more than one complex value w corresponding to the 
oscillation of the shell together with the adjacent fluid layers. The remaining Eigen values found were purely imaginary. They 
correspond to periodic movements of a fluid with an almost motionless shell. 
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In @ 


Figure 2 Addiction lm @ from wave number k 


3. NUMERICAL RESULTS 


Consider the option of natural oscillations, when the shell is filled with liquid. Fig. 1 and 2 shows, respectively, the dispersion curves 


in the dependence of Re w, Im w on the wave number of the k-first mode, whose damping coefficients are the smallest, and their 


Eigen values can be complex. In accordance with the numbering of the graphs, four different values of the coefficient n were set: 1) 


0.0009: 2) 0.0018 3) 0.15 4) 0.018, with the remaining parameters according to (1). It follows that there is a minimum critical value of 


the viscosity coefficient nx, above which, in the zone of high wave numbers of the first mode, appear as an a periodic wave number. 


As a result of a numerical experiment, it was found that the critical values of the viscosity coefficient nx are in the interval 


[0.0120 0.0125]. 
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